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1.  INTRODUCTION 

v,In  previous  work,!-!,-?]  we  have  used  a  Nanopar 
resistive  pulse  analycer  to  determine  the 
temporal  evolution  of  the  distribution  of 
clusters  formed  by  antigen  coated  latex 
microspheres  cross  linked  by  the  complementary 
antibody.  This  antigen-antibody  system  is 
conceptually  isomorphic  to  an  organic  polymer 
reaction  involving  RAf  units  (latex  microspheres 
with  f  active  antigens)  cross  linked  by  Bt 
units  (complementary  antibody)  with  only 
A-B  bonding  allowed. 

The  form  of  the  cluster  size  distribution 
was  studied  as  a  function  of  the  time,  concen¬ 
tration  of  RAt-  units,  and  concentration  of 
B  j  units.  The  experimental  results  were 
interpreted  theoretically  in  terms  of  solu¬ 
tions  of  the  unidirectional  Smoluchowski 
equations.yphese  equations  are 

—  ■  2  n=l  n'n  •"  Xn-n'V 

(1) 

ann'  xnxn' • 

Here,  .Xn  is  the  mole  fraction  of  clusters 
containing  n  monomers,  and  the  annt's  are 
the  bimolecular  reaction  rate  coefficients. 
These  coefficients  contain  the  information 
which  fully  determine  detailed  kinetic  evolu¬ 
tion  of  the  cluster  size  distribution.  To 
be  specific  we  have  found  that  it  is 
theoretically  convenient  and  experimentally 
acceptable  to  regard  the  structure  of  the 
matrix  element  anni  as  separable  into  two 
factors  viz 

ann '  *  8  Ffn ,n ' )  •  (-) 

Here  F(n,n')  describes  the  dependence  of  a 
on  the  numbers  n,  and  n'  of  units  in  each 
pair  of  interacting  clusters.  F  is  dimension¬ 
less,  and  can  be  normalized  for  convenience 
by  the  condition 

Ffl.lW.  Ida) 

We  shall  denote  F(n,n')  as  the  cluster  reac¬ 
tion  factor.  The  factor  g  is  independent 
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of  the  size  of  the  interacting  clusters  and 
is  a  measure  of  the  "stickiness”  associ¬ 
ated  with  the  binding  of  one  single  unit 
(in  a  cluster  of  n)  to  a  single  complementary 
monomer  (in  a  cluster  n').  The  binding  factor 
g  is  a  reflection  of  the  specific  properties 
of  the  antigen-antibody  bond,  and  the 
conditions  of  the  solution  and  the  coated 
spheres  which  can  be  experimentally  varied. 

These  conditions  include  the  initial  number 
(f)  of  antigenic  sites  on  each  sphere,  the 
antibody  concentration  [ab0]  ,  the  concen¬ 
tration  of  microspheres  initially  placed 
in  solution  [CG]  and  the  pH  and  salt  concen¬ 
tration  of  the  reaction  mixture.  Thus  the 
bimolecular  reaction  rate  coefficient  anni 
is  a  rich  source  of  information  on  the  basic 
physical  factors  which  control  the  process 
of  aggregation.  The  factor  F(n.n')  provides 
information  on  how  the  size  of  clusters  at 
the  present  time  results  from  the  aggregat¬ 
ing  smaller  clusters  at  earlier  times. 
Furthermore,  as  we  shall  see,  F(n,n')  contains 
information  on  the  Hausdorff - dimensionabi  1  ity 
(D)  of  the  evolving  clusters.  Here  we  can 
regard  the  clusters  as  fractals,  and  the 
factor  F  is  a  measure  of  the  topological 
and  dynamical  factors  which  describe  the 
coalescence  of  developed  clusters. 

The  factor  g  on  the  other  hand  is  of  great 
practical  importance.  For  in  the  use  of 
latex  immunoassays  [3, 4, 5, 6]  it  provides 
a  quantitative  means  for  optimizing  the 
reagents  so  as  to  maximize  the  sensitivity 
and  the  speed  of  the  agglutination  reaction.  ’ 
It  also  enables  the  immunologist  to  determine  * 
such  useful  parameters  as  the  binding  constant 
[7]  for  the  interaction  between  antibody 
in  solution,  and  the  antigen  bound  to  the 
microsphere.  It  also  permits  a  determination 
of  the  antibody  lifetime  [2],  and  the  degree  ,,r 
of  coverage  of  the  sphere  by  antigen  |2,“]. 

Our  prior  work  [1,2]  established  the  form  i 
of  the  cluster  size  distributions.  Within 
the  experimental  error  the  distribution  was 
consistent  with  the  form: 

rxn/X,  ■  n-bl^nbin-i^  V^ 

Here  X  is  the  total  mole  fraction  of  monomer  J _ 
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made  tor  the  value  ot'  the  exponential 
The  best  values  of  found  in  each  ease 
were  as  follows: 


initially  added  to  the  solution,  and 

b  =  (1  -  j_Ln  )  (4a) 

X 

is  a  measure  of  the  degree  of  reaction.  For 
our  system 


0  \  b  <  1.  (4b) 

We  also  showed  tha.t  asymptotically  as  b-1 

( Xn/X 1  takes  the  form 

1  1  (5a) 

n  *  b  nT 

with  -*3/2  (5b) 

Finally  by  investigating  the  form  of  known 
analytic  solutions  of  the  Smoluchowski 
equation  we  reported  that  the  data  found 
experimentally  could  be  produced  if  the 

factor  F  had  the  form 

Ffn.n' )  =  — —  .  (6) 

In  the  present  communication  we  report 
that  we  have  obtained  [7]  computer  solutions 
of  the  Smoluchowski  equations.  These 
solutions  permit  us  to  investigate  the  effects 
of  a  much  broader  class  of  forms  for  the 
cluster  reaction  factor  F(n,n')  upon  the 
distribution  of  cluster  sites.  Also  we 
have  shown  that  physical  consideration  of 
the  processes  of  coalescence  of  clusters 
enables  us  to  make  a  theoretical  connection 
between  the  form  of  F  and  the  topological 
characteritation  of  the  clusters.  This 
argument  permits  a  deduction  of  the  Hausdorff 
dimensionality  of  these  clusters  of  -latex 
microspheres  cross  linked  by  antibody. 

:.  DETERMINATION  OF  POSSIBLE  FORMS  OF  CLUSTER 
REACTION  FACTOR:  CONNECTION  WITH  HAUSDORFF 
DIMENSIONALITY 

Previous  analytic  solutions  of  the 
unidirectional  Smoluchowski  equations  have 
been  obtained  using  the  following  forms 
for  ann i  (7,8) 


ann 1  *  A 

=  B(n*n') 

(a) 

fb) 

(7) 

=  Cnn ' 

(c) 

We 

have  broadened 

the 

class  of  such 

forms 

by 

considering  as 

well 

the  following 

three 

functions 

3nn'  1  g(nw»n'w) 

fa) 

or 

ann '  *  g)W--n"- 

/Zy 

fb) 

(8) 

or 

ann  '  1  g<  nn' 

(c) 

Using  the  procedures  described  briefly  below 
we  found  that  any  of  these  choices  for  anni 
could  fit  the  data  withtn  experimental  error 
provided  that  an  appropriate  choice  was 


Eqn  8(a) 
Eqn  8(b) 
iqn  8(c) 


a  *  1.05  t  0.04 
u  *  1.02  t  0.05 
a  =  0.99  t  0.05 


19) 


It  is  interesting 
these  three  forms 
different  when  n^n' 
terms  depend  upon 
viz 


to  observe  that  each  of 
for  the  ann.'s  are  quite 
however  their  diagonal 
n  in  an  identical  manner 


ann-v  n  J. 


(10) 


It  is  possible  to  obtain  a  physical  inter¬ 
pretation  of  the  form  of  Eqn  8c  in  the  follow¬ 
ing  manner  .  Let  us  regard  each  cluster 
as  a  fractal  whose  Hausdorff  dimensionality 
is  D.  Thus,  the  relationship  between  number 
n  of  units  in  a  cluster  and  the  effective 
radius  R  of  the  cluster  is  given  by 

n  -v  R°.  (11) 

Let  us  now  assume  that  two  clusters,  each 
having  n  and  n’  units  respectively,  can 
only  interact  and  bind  with  one  another 
by  cross  linking  of  those  units  near  the 
surface  of  each  cluster.  Suppose  that  those 
units  on  the  surface  are  contained  within 
a  thickness  AR  around  the  radius  R  and  that 
AR  is  independent  of  R.  Then  the  number 
of  units  (An)  in  the  cluster  of  n  available 
for  bonding  is 

An  -v  AR  I  12a) 

or  An  n,  RD-1&R  ( 12b) 


or  using  Eqn  11  we  can  write 

Ann.  nfD-n/0  4r. 


(15) 


Similarly  the  number  available  for  bonding 
on  the  n’  cluster  is 


An' 


'D-»/DAR. 


Since  anni  is  expected  to  be  proportional 
to  AnAn'  it  follows  at  once  that 


ann '  ^  (nn‘ 


.D-ll/D 


(  14) 


On  comparing  this  with  Eqn  8c  we  find 


•a  D-l 
2  *  D 


115) 


This  then  provides  a  relationship  between 
the  experimental  parameter  and  the  Hausdorff 
dimensionality  D. 

Using  ,)  «  0.99  t  0.05  we  conclude  that 
for  our  system  that 


i'  -  t  ;  (  lb) 

Theoretical  calculations  based  on  equilibrium 
considerations  by  Parisi  and  Sourlas  [10| 
and  F.  Family  (11)  have  shown  that  in  the 
case  of  lattice  animals  the  scaling  exponent 
v  defined  by  the  Eqn  R  v  nv  has  the  value 
1/d.  Since  D=l/v  for  lattice  animals  our 
present  finding  is  in  agreement  with  their 
result. 

The  argument  above  is,  of  course,  a  crude 
one  in  that  it  regards  a  cluster  of  n 
particles  as  having  an  active  region  AR 
which  is  independent  of  both  the  sice  of 

that  cluster  (nl  and  of  the  sice  of  the 
cluster  (n'l  with  which  it  interacts.  In 
general  this  may  not  be  the  case.  In  our 
vifew,  however,  it  seems  likely  that  at  least 
for  large  clusters  the  active  region  for 
cluster-cluster  penetration  may  indeed  satisfy 
our  assumption.  The  work  of  M.  Plischke 
and  2.  Race  reported  that  this  conference 
is  relevant  to  this  issue  in  the  case  of 

the  Eden  model  and  DLA  in  two  dimensions. 
Clearly  it  would  be  most  desirable  to 
investigate  theoretically  the  size  iR  of 
the  interaction  cone  for  large  clusters 
using  computer  simulation  for 

three-dimensional  model  systems. 

Finally  we  discuss  the  manner  in  which 
we  deduced  the  best  value  of  *  for  each 

of  the  three  choices  for  anni  as  given  in 

Sa,b,c.  The  determination  of  based  on 

the  fact  that  the  experimentally  measurable 
quantity  b(6)  as  defined  bv , 

£  Ci 

bf  61  =  1  —f, - - 

iici 

remains  below  a  maximum  value  which  is  a 
function  of  u  .  The  behavior  of  b(6)  as 

a  function  of  a  dimensionless  time  T  is 
shown  in  Figure  1  from  which  it  can  be  seen 
that,  for  each  „ ,  b{0)  does  not  monotonica 1 ly 
increase  as  a  function  of  T  but  instead 
reaches  a  distinct  maximum  value.  This 
maximum  value  of  bf6)  is  a  well  behaved 
function  of  „  as  is  shown  in  Figure  2.  This 
figure  shows  b(6)max  as  a  function  of  u  for 
each  of  the  functional  forms  of  the 
coefficients  for  which  numerical  solutions 
of  the  Smoluchowski  equations  were  obtained. 
This  functional  relationship  between  „  and 


substantial  traction  ot  measurements  i' 
large  values  of  I,  hence  this  drop  const  i  tut 
an  experimental  measurement  of  bi6)nax: 

bWmax  s  *4y  -  -01  • 

Figure  4  is  a  detailed  view  of  Figure 
Z  was  used  to  convert  this  measurement  bf6)-jax 
into  measurements  of  <«. : 

u >  =  1.05  ♦  .04  for 

at  =  1.02  z  .05  for  ti*/-.j-/2)2 

wi  »  .99  z  .05  for  (ij)w/2. 
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experimental  measurement  of  b(blnax. 

The  experimental  value  of  b(6)max  was 
determined  from  roughly  800  measurements 
made  of  the  cluster  size  distributions  under 
a  wide  variety  of  experimental  conditions. 

It  was  found  from  a  frequency  histogram 
of  the  number  of  measurements  made  of  given 
values  of  bib),  shown  here  as  Figure  3, 
that  the  number  of  times  values  of  b(6)  were 
observed  drops  dramatically  as  bib)  increases 
hevond  .49.  This  Jatn  set  included  a 


FIGURE  I 

The  maximum  values  of  b ( 6 }  are  shown  through 
graphs  of  b(6)  versus  T  from  kinetic 
coefficients  of  the  form  iu  ♦  (top)  and 
i }*  (bottom)  with  values  of  u  ranging  from 
0.0  (highest  curve)  to  1.0  (lowest  curve) 
in  increments  of  0.2.  For  each  form  of  the 
coefficients  and  each  value  of  w  these  curves 
exhibit  a  maximum,  above  which  no  value  of 
b(6)  should  be  observed. 


FIGURE  2 

The  similarity  for  asymptotic  shape  of  the 
distributions  as  predicted  by  three  different 
forms  of  the  kinetic  coefficients  is  shown 
by  a  graph  of  b(6)max  versus  uij.  The  three 
forms  of  the  coefficients  are  £•  ♦  j"  (solid 
line),  (ijW2  (dashed  line),  and  ( iw/2*  jw/2)  2 
(dotted  line). 


